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Based on a recently introduced operator algebra for the description of a class of 
CN ' integrable quantum liquids we define the ground states for all canonical ensembles of 

these systems. We consider the particular case of the Hubbard chain in a magnetic 
q^ ' field and chemical potential. The ground states of all canonical ensembles of the 

o" 

model can be generated by acting onto the electron vacuum (densities n < 1) or 
hole vacuum (densities n > 1), suitable pseudoparticle creation operators. We also 
evaluate the energy gaps of the non- lowest- weight states (non - LWS's) and non- 
highest- weight states (non - HWS's) of the eta-spin and spin algebras relative to 
the corresponding ground states. For all sectors of parameter space and symmetries 
the exact ground state of the many-electron problem is in the pseudoparticle basis 
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the non-interacting pseudoparticle ground state. This plays a central role in the 



ctf 



pseudoparticle perturbation theory. 
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I. INTRODUCTION 

In contrast to higher dimension, the many-electron one-dimensional problem is non per- 
turbative. However, it also represents one of the few cases of many-body theory where there 
exist a number of exact solutions available. The most important class of such electronic 
systems are the integrable models solvable by the "Bethe ansatz" (BA) [|I],0,[M|- Within 
these integrable "Luttinger-liquids" 0^] the multicomponent models describing interacting 
fermions have been extensively investigated PiPinflPI CTCTJ^l ■ in spite of the non- Fermi 
liquid behavior, these systems show in the sectors of parameter space of symmetry U(l) a 
"Landau-liquid" character gJT|Jll|Jl^Jl|Jl|Jl5 . 



For some particular choices of canonical ensembles of the Hubbard chain || and other 
integrable models it has been assumed that the ground state corresponds to filling sym- 
metrically around the origin the BA quantum numbers [I^ |l0||13i|l4||l6i|l7||18| . Following this 
assumption, the pseudoparticle-operator algebra introduced in Refs. |IT| , [13| , [r^| revealed that, 
in contrast to the electronic basis, the pseudoparticle basis is perturbative. In particular, 
it was found that the natural reference state of the pseudoparticle perturbation theory is 
the many-electron ground state, which in the new basis is the non-interacting pseudopar- 
ticle ground state. This implies that this state plays a crucial role in the pseudoparticle 
perturbation theory. 

In this paper we present a detailed study of the ground-state problem which confirms 
that the many-electron ground state corresponds to filling symmetrically around the origin 
the BA quantum numbers and that in the pseudoparticle basis it is the non-interacting 
pseudoparticle ground state. One of our goals is to deepen and generalize to all sectors 
of parameter space of the Hubbard chain in the presence of a magnetic field and chemical 
potential the study of that ground state. Our study refers directly to the canonical ensembles 
belonging to the four sectors of parameter space where the symmetry is U(l) ® U(l) which 
correspond to electronic densities < n < 1 and 1 < n < 2 and spin densities — n < m < 
and — < m < n (for < n < 1) and — (2 — n) < m < and < m < (2 — n) (for 



1 < n < 2). However, it also provides the results for the sectors of higher symmetry and, 
therefore, for all canonical ensembles of the model which correspond to electronic densities 

< n < 2 and spin densities — n < m < n (for < n < 1) and — (2 — n) < m < (2 — n) (for 

1 < n < 2). 

Previous studies of the Hubbard chain have focused mainly on electronic densities < 
n < 1 and spin densities < m < n. However, in spite of the fact that the symmetries of the 
model provide information for the values of the physical quantities at densities 1 < n < 2 
and spin densities — (2 — n) < m < (2 — n) and densities < n < 1 and spin densities 
— n < m < from the corresponding values at densities < n < 1 and spin densities 
< m < n, we show in this paper that from the point of view of the pseudoparticle 
operator algebra it is useful to consider explicitly all parameter-space sectors. 

In the case of integrable models of simple Abelian U(l) symmetry the elementary ex- 
citations are generated by a single type of pseudoparticles |l!|. In this paper we find that 
the description of the low-energy excitations of non-Abelian integrable Hamiltonians in- 
volves a larger set of pseudoparticles. We show that the description of all gapless excitations 
branches of the four U(l) ® U(l) sectors of parameter space of the Hubbard chain in a mag- 
netic field and chemical potential involves eight types of pseudoparticles. The corresponding 
pseudoparticle algebra generates the low-energy Hamiltonian eigenstates from the electron 
vacuum (densities < n < 1) and hole vacuum (densities 1 < n < 2), respectively. The 
related study of the low-energy physics of the sectors of parameter space of higher symmetry 
in terms of the above pseudoparticles will be presented elsewhere [p0| . This study leads to 



the symmetry transformations of the eight types of pseudoparticles and reveals that the 
holons, antiholons, and two types of spinons of Ref. |21| are nothing but limiting cases of 
our general pseudoparticles. 

In the pseudoparticle basis the ground state of the many-body problem is, for all canoni- 
cal ensembles and corresponding symmetries, a "non-interacting" state, i.e., a simple Slater 
determinant of filled pseudoparticle levels. In addition, in the above sectors of symme- 
try U{1) <S> U(l) the elementary excitations simply correspond to pseudoparticle-pseudohole 



processes around this ground state. 

The present study of the low-energy physics of the four sectors of symmetry U(l)<g>U(l) 
of the Hubbard chain in a magnetic field and chemical potential introduces the eight types of 
pseudoparticles needed for a deeper understanding and description of the gapless excitations 
in the sectors of higher symmetry 5*0(4), SU(2) ® U(l), and U(l) ® SU(2). Therefore, the 
present study is a necessary step for the full characterization and understanding of the 
low-energy physics of the integrable non-Abelian many-electron quantum problems. 

In Section II we generalize the perturbation theory introduced in Ref. |TTi|l3| to all the 
U(l) ® U(l) sectors of parameter space of the model. This requires the introduction of eight 
pseudoparticle branches. 

In Sec. Ill we confirm that in the case of canonical ensembles belonging the sectors of 
parameter space of lowest symmetry U(l) <S> U(l) the ground state corresponds to filling 
symmetrically around the origin the BA quantum numbers. 

In Section IV we show that in all canonical ensembles of the four U(l) ® U(l) sectors of 
parameter space the non-LWS's and non-HWS's of the eta-spin and spin algebras have an 
energy gap relative to the corresponding ground state. In addition, we evaluate the smallest 
of these gaps. 

Finally, in Sec. V we present the discussion and concluding remarks. 

II. PSEUDOPARTICLE BRANCHES OF THE FOUR 17(1) ® 17(1) SECTORS 
Consider the Hamiltonian describing the Hubbard chain in a magnetic field H and chem- 
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H = H SO (4) + 2\xr\ z + 2^HS Z , (1) 

where 

Hsow = ~tJ2 [clc J+la + c] +la c ja ] + £/£[c] T c,. T - l/2][c] [Cji - 1/2] . (2) 



3,0 



Here the operator c] CT (cj ct ) creates (annihilates) one spin a electron at the site j and t, U, fi, 
H, and Hq are the transfer integral, the onsite Coulomb interaction, the chemical potential, 
the magnetic field, and the Bohr magneton, respectively. The operators, 

& = "[#.-£&], ^=~E^> ( 3 ) 



are the diagonal generators of the SU(2) eta-spin and spin algebras, respectively |22],|23 
a refers to the spin projections a =f , J, when used as an operator or function index and 
is given by a = ±1 otherwise. In Eq. (3) N ff = Y.j Cj a Cj a is the number operator for er 
spin-projection electrons. 

The model (1) — (2) describes an interacting quantum system of iVj up-spin electrons 
and iVj_ down-spin electrons on a chain of N a sites with lattice constant a. We use periodic 
boundary conditions and consider N a to be even and, when N = N a (half filling), both iVj 
and Ni to be odd. Henceforth we employ units such that a = t = fio = h= 1. The Fermi 
momenta are given by kp a = irn a and kp = [&ft + &fj.]/2 = irn/2, where n a = N a /N a and 
n = N/N a . The dimensionless onsite interaction is u = U /At. 

In the absence of the chemical-potential and magnetic- field terms the Hamiltonian (1) 



reduces to (2) and has SO (A) = SU(2) <8> SU(2)/Z 2 symmetry ||,|3|,||,|5],||,|3 . Since 
N a is even, the operator fj z + S z (see Eq. (3)) has only integer eigenvalues and all half-odd 
integer representations of SU{2) ®SU{2) are projected out |H|,^J. The two 577(2) algebras 
- eta spin and spin - have diagonal generators given by Eq. (3) and off-diagonal generators 
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V = E(-l)^T9l , V' = E(-l) J sV]t . ( 4 ) 

3 3 

and 

s = E c h c 3i » & = E c )i c 3i » (s) 

3 3 

respectively 

In the presence of both the magnetic field and chemical potential terms, the symmetry 
is reduced to U(l) ® U(l), with fj z and S z commuting with H. The eigenvalues rj z and S z 



are determined by the values of the conserved numbers, as shown by Eq. (3). According to 
these eigenvalues, the system has different symmetries as follows fT0| , [T3"|| : when r] z ^ and 
S z ^ the symmetry is £7(1) ® £7(1), for r] z = and S z ^0 (and \t = 0) it is 577(2) ® 17(1), 
when Tj z 7^ and 5 2 = it is £7(1) <8> 577(2), and at ?? 2 = and 5 2 = (and fi = 0) the 
symmetry is £0(4). The £7(1) (g> 77(1) symmetry sectors always correspond to two non-zero 
eigenvalues of the diagonal generators, whereas in sectors of higher symmetry, one or both 
of these eigenvalues vanish. There are four 77(1) <g> 77(1) sectors of parameter space which 
correspond to rj z < and S z < 0; t] z < and S z > 0; r] z > and S z < 0; and r\ z > and 
5 Z > 0. There are two 5*77(2) <g> 77(1) [and 77(1) eg) 577(2)] sectors of parameter space which 
correspond to S z < and 5 2 > (and to r\ z < and r^ 2 > 0). There is one 50(4) sector of 
parameter space [which is the above t] z = (and \x = 0) and S z = "point"]. 

We restrict our analysis to the four sectors of lowest symmetry 77(1) <g> 77(1). In this 
case we can define the parameters / = sgn(r] z )l and /' = sgn(S z )l which classify these four 
sectors: we denote each of them by (/,/') sector. The (—1,-1); (—1,1); (1,-1); and (1,1) 
sectors refer to electronic densities and spin densities < n < 1 and 0<m<n;0<n<l 
and — n < m < 0; 1 < n < 2 and < m < (2 — n); and 1 < n < 2 and —(2 — n) < m < 0, 
respectively. 

In the (—1,-1) sector the BA solution refers only to the LWS's of both the eta-spin and 
spin algebras. Note that the Hamiltonian eigenstates of the remaining three (I, I') sectors 
can be generated by multiple application of the operators tf (4) and 5"^ (5) to the LWS's of 
the (—1,-1) sector PJ2~3"|[ . We find below that in the sectors where r] z > and (or ) S z > 
the ground states and low-energy Hamiltonian eigenstates are HWS's of the eta-spin and (or 
) spin algebras. However, the generation of these HWS's from the corresponding LWS's with 
the same rj and 5 values [and r\ z = —r\ z and (or ) S z = —S z ] represents a too complicated 
and indirect description of the low-energy physics of the rj z > and (or ) S z > sectors of 
parameter space. Instead, in this paper we describe the low-energy physics of these sectors 
directly from the corresponding BA solutions. We emphasize that the set of LWS's and (or 
) HWS's of each of the four (I, V) sectors of symmetry 77(1) <g> 77(1) represent four alternative 
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choices for starting states to construct a complete set of Hamiltonian eigenstates [see Ref. 
|23| where the LWS's of the (—1,-1) sector are used as starting states]. 

Since in the (1,-1) [and (—1, 1)] sector the low-energy physics is determined by HWS's 
[and LWS's] of the eta-spin algebra and LWS's [and HWS's] of the spin algebra, hence- 
forth we refer the LWS's and (or ) HWS's of the four sectors as [LWS,LWS]'s [sector 
(-1,-1)], [LWS,HWS]'s [sector (-1,1)], [HWS,LWS]'s [sector (1,-1)], and [HWS,HWS]'s 
[sector (1, 1)]. (This notation refers to the order [eta-spin, spin].) 

Let 1 77s, S z ) be a regular BA Hamiltonian eigenstate corresponding to the canonical en- 
semble of eigenvalues rj z and S z . By regular BA Hamiltonian eigenstates we mean here these 
states to which the BA solution of any of the four U(l) <8> U(l) sectors directly refers. Such 
states can be [LWS,LWS]'s, [LWS,HWS]'s, [HWS,LWS]'s, or [HWS,HWS]'s according to the 
sector of parameter space. The BA solution for the sectors (—1,1), (1,-1), and (1,1) is 
obtained as for the usual sector (—1,-1), except that instead of the equations 

v\Vz,S z ) = S\rj z ,S z ) = 0, for r, z < , S z < , (6) 

which refer to [LWS,LWS]'s, the following equations have to be fulfilled 

ri\7 lz} S z } = &\r ]z ,S z } = 1 for Vz <0,S z >0, (7) 

for the (—1, 1) sector where the regular BA Hamiltonian eigenstates are [LWS,HWS]'s; 

fjHvz,S z ) = S\ Vz ,S z ) = 0, for Vz >0,S z <0, (8) 

for the (1, —1) sector where the regular BA Hamiltonian eigenstates are [HWS,LWS]'s; and 

V [ \Vz,S z ) = &\7] z ,S z ) = i for Vz >0,S z >0, (9) 

for the (1, 1) sector where the regular BA Hamiltonian eigenstates are [HWS,HWS]'s. 



In References [[TT],[U||T4]] the LWS's of the usual (—1,-1) sector were devided into two 
types: the LWS's I, which refer to real BA rapidities; and the LWS's II. Some or all the 



rapidities which describe the latter states are complex and non real fll3| . In the canonical 



ensembles of that U(1)®U(1) sector both the non-LWS's and the LWS's II have finite energy 
gaps relative to the ground state flcSfl . 

As in the case of the [LWS,LWS]'s of the (-1,-1) sector, the [LWS,HWS]'s, 
[HWS,LWS]'s, and [HWS,HWS]'s can be devided into [LWS,HWS]'s I, [HWS,LWS]'s I, and 
[HWS,HWS]'s I, described by real rapidities, and [LWS,HWS]'s II, [HWS,LWS]'s II, and 
[HWS,HWS]'s II. Some or all the rapidites which describe the latter states are complex and 
non real. In the four (/, /') sectors the non-LWS's and non-HWS's and the states II have 
an energy gap relative to the ground state, which is always a [LWS,LWS] I, [LWS,HWS] 
I, [HWS,LWS] I, or [HWS,HWS] I. In Sec. IV we will show that for the non-LWS's and 



non-HWS's, whereas the gaps of the states II will be evaluated in a following paper p8[ . At 
energy ranges smaller than those gaps the Hilbert subspace spanned by the regular BA states 
I with common rj z and S z eigenvalues coincides with the full Hilbert space of the quantum 
problem (spanned by states with these eigenvalues). The perturbation theory introduced 
in Refs. flTT| , ri"3| for the (—1, —1) sector is generalized here for all (1,1') sectors and refers to 
that Hilbert space. It corresponds to energy scales smaller than the above gaps. 

Henceforth, we will call state I any [LWS,LWS] I, [LWS,HWS] I, [HWS,LWS] I, 
or [HWS,HWS] I, and state II any [LWS,LWS] II, [LWS,HWS] II, [HWS,LWS] II, or 
[HWS,HWS] II. We note that a state I or a state II is always a regular BA Hamiltonian 
eigenstate and, therefore, is not, by definition, a non-LWS and (or ) non-HWS. 

Let us generalize to the states I of all the U(1)®U(1) sectors the pseudoparticle-operator 
algebra of the (—1, —1) sector. In each of the (I, /') sectors we have two types of pseudopar- 
ticles. Combining the four sectors we have eight types of operators b * „ j« and b qa (i,i') which 
obey the usual fermionic algebra 

{bqa(l,l')i bq' a '(l",l">)} = ^q,q'^a,a'^l,l"^l',l'" , (10) 

and 

i b qa(l,l>)i b q'a'(l",l'")} = ° > V>qa(l,l'), \'<*' Q" ,*'") } = ° ' ( ll ) 



Therefore, there are eight quantum numbers or "colors" a(l,l'), which are c(— 1, — 1), 
c(l, — 1), c(— 1,1), c(l,l), s(— 1, — 1), s(l, — 1), s(— 1,1), and s(l,l). These correspond to 
eight different branches of a(l, I') pseudoparticles. 

In each (I, I') U(1)®U(1) sector only the corresponding a(l, I') pseudoparticles participate 
in the construction of the states I, as we discuss below. Therefore, the algebra of Eqs. 
(10) and (11) may be replaced by {^^^.Va'fM')} = <W<W and {b\ a{lV yb\, al{lV) } = 0, 
respectively. The corresponding discrete pseudomomentum values are 



'a 



r a(l,l') 



where Ij are consecutive integers or half integers. The representation (12) was used by 
Yang and Yang in the case of the one-dimensional boson gas with repulsive delta-function 



interaction 29 



There are N*^^ possible I" ' values. A state I is specified by the distribution of N a (ij>) 
occupied values, which we call a(l, I') pseudoparticles, over the N*, t j« available values. 

Since only single and zero occupancy of the values Ij ' are allowed, only pseudoparticles 
of the color a(l,V) can occupy the states labeled by the numbers I" . Therefore, the 
pseudoparticles have a fermionic character, as assured by the anticommuting algebra (10) 
and (11). (The BA wave function vanishes for configurations showing other than single and 
zero occupancy of the values /• ' |)4|,|l3|.) For the Hubbard model the BA spatial wave 
function for the states I depends on the quantum numbers Ij ' through two sets of real 
numbers, which many authors call rapidities. 

The BA solution for the model (1) can be obtained directly for all sectors of parameter 
space. In Appendix A we provide a short presentation of that solution in the case of the 
present four (1,1') sectors of symmetry U(l) <8> U(l). In that Appendix we introduce the 
rapidity functions, two-pseudoparticle phase shifts, and some other quantities needed in 
this paper. Equations (Al) and (A2) define the above rapidites for a given choice of the 
distribution occupancy of the numbers /■ ' , i.e. for a given state I. (For the states II some 
or all rapidities are complex, non-real, numbers |[28|| .) 



There are N*,^ — N a Q^ empty values, which we call a(l, I') pseudoholes. Carrying the 
BA solution in the same manner as for the usual (—1,-1) sector we find the numbers N*, t ^ 
and N a ny\ for the remaining three sectors. In Table 1 we give these numbers for the four 
sectors. These are conserving numbers involving the numbers of lattice sites N a and of o 
electrons N a . 

The eta spin rj and spin S values [and the corresponding eigenvalues of the diagonal 
generators (3), r\ z = If] and S z = I'S] can be expressed in terms of the numbers of pseudoholes 
as follows 

V = \\K m ~ N cm ] , (13) 

and 

S = \[K m - N. m ] = \[N cm - 2N sm ] , (14) 

respectively. The numbers Ij ' are integers (or half integers) for iV^/) even (or odd), and 
Ij ' are integers (or half integers) for N*^^ odd (or even). 

Let \V; — 1, =Fl) and \V; 1, =Fl) be the vacuum states which correspond to taking the limit 
of the electronic density n — > and n — > 2, respectively, keeping ±iV| > ±N±. In \V; —1, =Fl) 
and \V; 1, =Fl) there are only holes and electrons, respectively. 

Since the colors a(l, V) and the pseudo momentum q are the only quantum numbers 
involved in the description of the pseudoparticles whose occupancy configurations define the 
states I, all the corresponding allowed configurations can be generated by applying to the 
corresponding vacuum \V;l,l') the J2a^a(i,i') creation operators b qa n v -, : i.e. 

qN a(l I') 

\V Z ,S Z )= I][ II bl m ]\V;l,l'), (15) 

a=c,s qj=qi 

[here / = sgn(r) z ), I' = sgn(S z )] where the set of qj values, qi, q 2 , ■■■, QN a(lll) , refer to the 
N a (i,v) occupied pseudomomenta out of the N*, t v \ available pseudomomentum values. 

In each canonical ensemble belonging the (/, /') sector, out of the total of regular BA 
Hamiltonian eigenstates I and II, 
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N a 

. N ki,n . 




N s (i,i') 


+ 



N a 
N s (i,i>) — 2 



N a 



N; m + 1 



N a 



N kun ~ l 



N a 



N, 



s(l,V) 



(16) 



(this formula is a generalization of the (—1,-1) result of Ref. ||23||), there are 



N„ 



N S {i,v) + N; m 



N* 



N. 



s(l,l') 



(17) 



Hamiltonian eigenstates I. (The square brackets in the above equation refer to the usual 
combinatoric coeffi cents.) 

We emphasize that in the suitable limits formula (17) also applies to canonical ensembles 
of symmetry SO(A), SU(2) <g)U(l), and U(l) ®SU(2). For instance, in the n = 1 and m = 
SO (4) canonical ensemble for any of the four choices of (/, /') numbers Eq. (17) gives one. In 
this case these four choices are alternative representations of the same canonical ensemble, 
i.e. when rj z = (and fi = 0) and S z = there is only one state I with rj = S = 0. This 
is both a LWS and a HWS of the eta-spin and spin algebras and is the SO (4) ground state 

. In this canonical ensemble there are neither LWS's I nor HWS's I excited singlet states 



of the eta-spin and (or ) spin algebras. 

In each canonical ensemble of eigenvalues r] z and S z the states I of the form (15) and 
of total number given by Eq. (17) constitute a complete orthonormal basis which spans an 
important Hilbert subspace, which we call Tij. At energy scales smaller than the gaps for 
non-LWS's, non-HWS's, and states II, Tij represents the full accessible Hilbert space. We 
emphasize that, in general, many states I in Tii have energies larger than these gaps. The 
low-energy physics is determined only by the states I whose energies are smaller than such 
gaps. 

In Ref. |Tj|, and following Refs. P Jl~0| , pT|rf7| , p!8[] , it was assumed for the case of the sector 
(—1, —1) that out of all corresponding states I of form (15), the ground state of eigenvalues rj z 
and S z corresponds to filling symmetrically around the origin N a (-i-i) consecutive I" 
values of all colors a(— 1, — 1). In the present general case this leads to 



ra(-l,-l) 
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q (+) l , 
\0;v z ,S z )= ni IT ^(M')]l y ;^'), (18) 

a=c,s (-) 

9=9Fa((,i') 

where when N a ^^ is odd (or even) and /" ' are integers (or half integers) the pseudo- Fermi 
points are symmetric and read 

v { F* m = -iF* { i,n = w [Na ^ ~ 1] ■ (19) 

If all pseudo-Fermi points are symmetric the state (18) has zero momentum and is nonde- 
generate. On the other hand, when at least one of the pseudoparticle or antipseudoparticle 
numbers N a (i^} is odd (or even) and Ij ' are half integers (or integers) the corresponding 
pseudo-Fermi points are nonsymmetric and read either 

<lF* m = TrN am , -gfcV) = ^- [N am - 2] , (20) 

or 

tilm = ^K(M') - 2] , -QFam = Tr^m ■ (21) 

In this case the state (18) has finite momentum and is degenerate. Equivalent expressions 

can be obtained for the limits of the pseudo-Brioullin zones, gL/j/), if we replace N a (i t i^ 

by N*njn in Eqs. (19) — (21). Except for terms of order l/N a , we have that q F Lw) — 
,(-) ,^ „(+) - „H 



^k(M') = ^W) and 9a(M') = -9a(/,J') = 9a(J,J'), where 



9Fa(J,I') - 77 . Qa(l,l') ~ 77 • (^J 

In some studies the full expressions (19) — (21) have to be used because the terms of order 
1/N a play an important role. On the other hand, many quantities are in the thermodynamic 
limit insensitive to these 1/N a corrections and we can replace q F ^n t i\ in many expressions 
by the pseudomomenta ±qp a (i,i') (22). In Table 1 we present the values of the pseudo- Fermi 
points and limits of the pseudo-Brillouin zones (22) for the four (I, I') sectors. 

In Sec. Ill we will confirm that the ground state of canonical ensembles of sectors of 
parameter space of symmetry U(l) <S> U(l) has the form (18). 
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The Hamiltonian eigenstates (15), in number of (17) and of which the state (18) repre- 
sents a particular choice, can be rewritten relatively to the latter state. In this case they 
correspond to pseudoparticle-pseudohole processes around the reference configuration (18) 
and read 

N a(_l,l') 

ph 

Ks M )= Hi II fcUwWdl ;^^}, (23) 

a=c,s ij=l 

where qj (<&) defines the different locations of the pseudoparticles (pseudoholes) relatively 
to (18) and iV°jJ ' is the number of a{l,V) pseudoparticle-pseudohole processes. When the 
starting ground state and the final state of an excitation I have different electronic numbers, 
N a , it can be decomposed in a ground-state - ground-state transition and a pseudoparticle- 
pseudohole excitation around the end ground state [|(J . 



All the Hamiltonian eigenstates (15) and (23) are states I. Non-LWS's and non-HWS's 
are generated by acting raising or lowering generators (4) and (5) on these states. 

In the Hilbert subspace Tii spanned by the states I the Hubbard model can be written in 
the pseudoparticle basis. The derivation is as for the (—1,-1) sector and the Hamiltoniamn 
(1) reads [|TJ 

H = H so{4) + |/i| £[1 - #c(mo(<z)] + Vo\H\ £[1 - N sm (q)} 
1 i 

= H SOW + |//| £[1 - N cm (q)} + ^o\H\[J2N cm (q) - 2^iV s(M0 (g)] , (24) 

<? ii 

where the Hamiltonian Hso(4) (2) is of the form 

Hso{A) = E [N cm (q){-2tcos[K^(q)} - U/2} + VIA] , (25) 

i 

Kiji(q) is a rapidity operator [|13|] whose eigenvalues, Kiy(q), are studied in Appendix A [see 

Eqs. (A3) and (A4)], and N c Q t i>)(q) is the c(l,l') pseudomomentum distribution operator. 

The operator N a (iji)(q) has the form 

N a (i,i>){q) = &Ja(J,i') 6 «°(',i') ' ( 26 ) 

and we can write the a{l,V) pseudoparticle number operator as 
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#*(!,«')= E^aWoGz)- (27) 

It follows that the operators rj z and S z read 

7), = /£[1 - ^c(M')(?)] , S* = l 'J2l l - N. m (q)] ■ (28) 

q q 

This, together with the following relations valid for all (I, V) sectors, 

-Ifi = |/i| ; -I'H = \H\ , (29) 

justifies the form (24) of the Hamiltonian (1) in the pseudoparticle basis. 

In normal order relatively to the state (18) the expression of this Hamiltonian involves 
the normal-ordered operators ||i~3| . |i4" 



: N qa{hV) :-: b qa{lv) b qoi{hV) : , (30) 

and has pseudoparticle forward-scattering terms only. (The expression of the normal-ordered 
operator : b „ v \b qa (i t i') : of the right-hand side (rhs) of Eq. (30) is given in Eq. (43) below.) 
As in the case of the (—1,-1) Hamiltonian of Ref. ||13| , it has an infinite number of terms 
which correspond to increasing scattering orders. To second order we find for each (/, V) 
sector 

: H :=E e «(M')W : N qam : +^EE faa'iliV') ■ N qa (i, V ) :: N q , a > m : + (31) 

a,q Z a,q a ',q' 

where the pseudoparticle bands e a (i,i>)(q) and the / functions are evaluated as in the case 
of the (—1,-1) sector |T3|JT8|| . Since the colors a(l,l') and a'(l,l') of two-pseudoparticle 
quantities, such as the above / functions and below phase shifts (see also Appendix A), 
refer to the same (I, V) numbers, in order to simplify our notation we call in this case 
these colors a and a', respectively. (The numbers (I, V) need to appear only once in these 
functions.) The expressions of the bands are 

e c (M')(<?) = CcCMofe) + [M - U M - ^1^1 ' ^ 32 ) 

and 
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e s (i,i'){q) = e° (M0 (g) + 2fi \H\ , (33) 

where the pseudoparticle bare spectra e° c n v s (q) and e% m (p) are given by 

e° cm (q) = -2tcoaK$(q) + 2t /° M ' dk*g (k,K$(q)) sink, (34) 



and 



e\ m (q) = 2t f Qu ' dkQg (k, S® (q)) sin k , (35) 

respectively. The phase shifts $^, a , are defined in Eqs. (A39)-(A42) and the functions K\jj{q) 
and S\i,{q) are the solutions of Eqs. (A3)-(A4) for the reference state (18). The parameters 
Qi : i> are given in Eq. (A22). The functions K\ii{q) and SiJ(q) can be defined in terms of 
the phase shifts (A39) and (A40), respectively, as follows 

K®{q) = q - f Qu ' dk*g (k, <>(,)) , q = f^' dh®£ (k, S$(q)) . (36) 

J -Qi,i< v J J -Qi,i> v ' 



The pseudoparticle bands e a (_i ! _i)(g) = e a (q) are plotted in Figs. 7 and 8 of Ref. |L7| 



In the present (I, I') sectors of parameter space of symmetry U{1) <8> U(l) the pseudoparticle 
energy spectra €c(i,i'){q) an d e s(i,i')(q) defined by Eqs. (32) and (33) vanish at the pseudo- 
Fermi-points, i.e., 

e «(i,i')(<l { F*(i,i')) = °- ( 37 ) 

The / functions, fa a *{q, q'), of the rhs of Eq. (31) read 

fZ'iq^q') = 27tv a ^i l) (q)^ l £ / (q,q')+27rv al{l>ll) (q')^ a (q',q) 

+ 12 12 27rv a , im $ l a i a (jq Fa ,, mj q)$ l a i a ,(jq Fa » {l j, )j q'), (38) 

j'=±l a"=c,s 

where the two-pseudoparticle phase shifts ^^ a ,(q,q') are defined by Eqs. (A18)-(A21) and 
the pseudoparticle group velocities are given by 

Va(i,i>){q) = -j 1 — • (39) 

The "light" velocities 
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v a (i,i>) = v a (i yV) (q Fa (i yV) ) , (40) 

play a determining role at the critical point and appear in the conformal-invariant expressions 



ll|JT3|| . (The velocities iW-i -1) = v a are plotted in Fig. 9 of Ref. [|17 |. 



Equation (31) is a generalization of the corresponding (—1, —1) Hamiltonian of Ref. [13 
for the three remaining sectors. We emphasize that in the present U(l) <g> U(l) sectors of 
parameter space and at energy scales smaller than the gaps for the non-LWS's, non-HWS's, 
and states II Eqs. (24), (25), and (31) refer to the expression of the full quantum-liquid 
Hamiltonian. In the electronic basis this is given by Eqs. (1) and (2). (In the case of the 
sector (— 1, —1), the pseudoparticle-operator representation (31) leads, in a natural way, to 
the low-energy spectrum studied in Refs. [P|,P|, pX>| , P~7| , [TB| ] . ) 



Both the / functions (38) and all the remaining higher order coefficients have universal 
forms in terms of the two-pseudoparticle phase shifts and pseudomomentum derivatives of 
the bands and their higher order "velocities". The two operators of the rhs of Eq. (31) are 



the Hamiltonian terms which are relevant at low energy [Tl]|T3[ . 

The perturbative character of the pseudoparticle basis follows from the fact that, in con- 
trast to the two-electron forward scattering amplitudes and vertices, the two-pseudoparticle 
/ functions (given by Eq. (38)) and the corresponding two-pseudoparticle forward-scattering 
amplitudes, which were calculated in Ref. [Kj for the (—1, —1) sector, do not diverge and 



are finite. 

The combination of Eqs. (32), (33), and (37) allows the derivation of the density and 
magnetization curves which read 

V = -*[-o ~ e °(z,z')(<?F C (M')) \x -2 e °(U')(^W))] \x , (41) 

and 

H = n 4 w (aw)) n (42) 

Zfi 
respectively, where, depending on the constraints imposed to the system, either x = H or 
x = m and either y = \x or y = n, respectively. 
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The excited states I, which in the canonical ensembles of the four sectors of lowest 
symmetry are the only gap less states, never involve pseudoparticles with different (I, I') 
numbers. In each of the four (I, I') sectors of symmetry U(l) <8> U(l) there is a selection rule 
which states that out of the eight pseudoparticle branches only pseudoparticle-pseudohole 
transitions in each of the two corresponding bands c(l, V) and s(l, I') are allowed. Therefore, 
in the canonical ensembles of the (I, V) sectors the states I (15), in number given by Eq. (17), 
refer to the two c(l, I') and s(l, I') bands only. This is also true for all states I associated 
with higher symmetry sectors ||20|| . 



III. THE GROUND STATE IN THE 17(1) ® 17(1) SECTORS 

In this section we use the pseudoparticle basis and the associated perturbative character 
of the quantum problem to justify that in canonical ensembles of the (I, I') sectors the ground 
state has the form (18) and corresponds to filling symmetrically around the origin the BA 
quantum numbers. 

To show that at given eigenvalues r] z and S z and within all the states I of form (15) and 
in total number given by Eq. (17) the Hamiltonian eigenstate (18) (which is degenerate 
when its momentum is finite) is the state of minimal energy, we use the pseudoparticle basis 
of Tii, where the Hamiltonian (1) and (2) is given by Eqs. (24) and (25), respectively. 

In Appendix B we study the general energy expression for all states I of Tii (with common 
r] z and S z eigenvalues), which is given by Eqs. (Bl), (B2), and (B4). In the case of zero-spin 
density, m = 0, the energy Egou) + (U/2)[N — N a /2] was studied and plotted by Shiba 
|T6| [here E° so ,^ is the energy (B2) and (B4) of the sector (—1, —1)]. In that Appendix we 
devote particular attention to the states I of minimal and maximal energies. The study of 
the energy expressions for the different Hamiltonian eigenstates with common eigenvalues 
7] z and S z identifies these two states I. The results are: 

First, the study of the energy (Bl) reveals that the energies of the states I of form (15) 
and with common eigenvalues r\ z and S z correspond to a continuous distribution without 
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energy gaps. 

Second, the minimal and maximal energies of that continuous distribution of energies of 
states with common eigenvalues rj z and S z corresponds to the states I of pseudomomentum 
distributions N a ^)(q) = N^^(q) and iV a(M /)(g) = iV* (M/) (g), respectively, where -/V° (M /)(<?) 
and N^n^Jq) are given by Eqs. (B3) and (B5), respectively. 

Third, since the energies of the states I of form (15) and with common eigenvalues i] z and 
S z correspond to a continuous distribution without energy gaps, in order to show that (Bll) 
and (B12) are the minimal and maximal values for these energies, respectively, it is enought 
to show that the energies [E — Eq] and [E — E*\ of the sub class of these states that can 
be generated from (18) and (BIO), respectively, by changing the distribution occupancies 
of an arbitrary small density of pseudoparticles relative to the distributions (B3) and (B5), 
are positive and negative, respectively. Therefore, within all states (15) in number of (17) 
it is enough to evaluate the energies of the Hamiltonian eigenstates which differ from the 
reference distributions (B3) and (B5) by changing the occupied pseudomomenta by a small 
density for the two branches of c(l, I') and s(l, I') pseudoparticles. 

Fourth, and in contrast to the electronic basis, we can use the perturbative character of 
the pseudoparticle basis pi to expand the energies [E — E ]/N a and [E — E*]/N a of the 



above sub class of states in the densities of excited «(/,/') pseudoparticles, n"( 1 ' 1 \ relative 
to the reference distributions (B3) and (B5). Furthermore, we can write the Hamiltonian 
(1) with (2) given by (25) in normal order relative to the reference states (18) and (BIO) 
|13|| . In Sec. II we have presented the normal-ordered Hamiltonian relative to the state (18), 



Eq. (31). The normal-order character of these Hamiltonians implies that the corresponding 
energies are given relative to the states (18) and (BIO), respectively. These Hamiltonians 
have an infinite number of terms [see Eq. (31) for the case of the state (18)] and, therefore, 
the evaluation of these energies seems to require the evaluation of an infinite number of en- 
ergy contributions. However, the perturbative character of the pseudoparticle basis makes 
the problem much easier. This perturbative character rests on the fact that the evaluation 
of these energies up to the i th order in these densities requires considering only the corre- 



sponding Hamiltonian terms of scattering orders less than or equal to i ||13|| . This follows 
from the linearity of the density of excited a(l, I') pseudoparticles, which are the elementary 
"particles" of the quantum liquid, in SN a ^ t i'\(q) = (r) z ,S z \ : N a (i,i')(q) : \i] z ,S z ). Here, the 
normal-ordered pseudomomentum distribution relative to the states (18) and (BIO) is given 
by 

: N am (q) := b\ a{lv) b qa{hV) - N° m {q) , (43) 

and 

: N am (q) ■= b\ am b qam - N* m (q) , (44) 

respectively. 

Fifth, the perturbative character of the pseudoparticle basis together with the above 
analysis implies that in each canonical ensemble it is enough to consider the states called, 
in the case of a (—1, —1) state (18), (B) and (C) in Refs. |H,I3]|. These states can also be 



defined in the general case of the (I, I') state (18) and of the state (BIO) and correspond to a 
small density of pseudoparticle-pseudohole processes relative to the reference distributions 
(B3) and (B5). (The states (A) of Refs. |TT| , |T3"|] are associated with changes, AN a , in the 



number of a electrons.) Moreover, one needs to evaluate the energies of these states up to 
second order in the density of a(l, I') pseudoparticles, n^ l ' l '\ only. This involves only the one- 
and two-pseudoparticle terms of the corresponding normal-ordered Hamiltonian. [For the 
case of the state (18) see Eq. (31).] We find that for the reference states (18) and (BIO) all 
such energies are positive and negative, respectively. (For the case of the (—1, —1) reference 
state (18) these energies were evaluated in Ref. fl3fl , and equal the energy spectrum studied 
in Ref. || when the number of both up-spin and down-spin electrons is kept constant.) 
This confirms that in each canonical ensemble of a (/,/') sector (18) and (BIO) are, within 
all states I of form (15), with common eigenvalues 7] z and S z , and in number of (17), the 
Hamiltonian eigenstates of minimal and maximal energies, respectively. This implies that 
the energies E, Eq. (Bl), of all remaining states I of that ensemble are such that 
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E <E <E*, (45) 

where Eq and E* are the energies (Bll) and (B12), respectively. In addition, for each canon- 
ical ensemble of eigenvalues r\ z and S z , [E* — E ] gives the energy width of the continuous 
distribution of energies corresponding to the whole set of (I, I') states I of form (15). 

We emphasize that only writing the Hamiltonian in the pseudoparticle-operator basis 
introduces the perturbative character of the quantum problem which simplified the above 
analysis. Here we have compared the energies of the states I only. Therefore, a complete 
proof of (18) being the Hamiltonian eigenstate of minimal energy in canonical ensembles 
of the present lowest-symmetry sectors of parameter space requires the evaluation of the 
energy gaps of the non-LWS's and non-HWS's and of the states II relative to that state, 
which have to be finite. While the gaps of the non-LWS's and non-HWS's are evaluated in 
Sec IV, the gaps of the states II are calculated elsewhere |2^ . 



Finally, the study of the spectrum of the states II of the sectors of higher symmetry 



SO (4), SU(2)®U(1), and U(1)®SU(2) p0p8j reveals that in canonical ensembles of these 



sectors the ground state is also a state I of form (18), as we discuss in Ref. ||20|| . [When 
rj z = and (or ) S z = 0, the a(±l,/') and (or ) a(l, ±1) pseudoparticles correspond to 
alternative representations of the same ground state.] Therefore, in all canonical ensembles 
of the quantum problem and in the pseudoparticle basis the ground state is a simple Slater 
determinant of pseudoparticle levels of the universal form given by Eq. (18). 

IV. ENERGY GAPS OF THE NON-LWS'S AND NON-HWS'S 

In Sec. Ill we have used the pseudoparticle basis to show that among all (/, V) states I 
with common eigenvalues r\ z and S z , in number of (17) and of the form (15), the Hamiltonian 
eigenstate (18) has minimal energy. Therefore, if in the corresponding canonical ensemble 
both all the non-LWS's and non-HWS's and all the states II (with common rj z and S z 
eigenvalues) have a gap relative to the state (18), this state is the ground state. In this 
section we show that all such non-LWS's and non-HWS's have an energy gap and calculate 
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the smallest of these gaps. 

Let us consider the Hamiltonian H so ^, Eq. (2). This represents the Hubbard chain (1) 
at zero magnetic field and chemical potential and commutes with the six generators of Eqs. 
(3) — (5) and, therefore, has 5*0(4) symmetry [^3|], as discussed in Sec. II. 



Let \t], S; ri z , S z ) be an arbitrary Hamiltonian eigenstate belonging the family of states 
with fixed eta spin and spin, r\ and S, respectively. In \rj, S; r/ z , S z ) r\ z and S z are the eta-spin 
and spin projections, respectively, of this particular state. The whole family of states with 
the same values of eta spin r\ and spin S but different eta-spin and spin projections can 
be generated from one of the four [LWS,LWS], [LWS,HWS], [HWS,LWS], and [HWS,HWS] 
with these values of eta spin and spin. The [LWS,LWS] starting state, for example, is the 
state \rj, S; — rj, — S) which belongs to the sector (—1,-1). The remaing three choices for 
starting states are \r], S; —r], S), \r), S;r], —S), and \rj, S;r], S), respectively. Depending on 
the LWS or HWS character of the starting state, the family of states is generated by acting 
operators f) and S or tf and ™ of Eqs. (4) and (5) onto that state. 

In Appendix C we consider fixed values of the chemical potential /z and magnetic field 
H such that /i ^ and H ^ 0. The signs of the chemical potential and magnetic field 
fix the signs of r\ z and S z and choose the particular (/,/') sector. According to Eq. (29), 
in the sectors (—1,-1), (—1,1), (1,-1), and (1,1) the chemical potential /i and magnetic 
field H are such that fi > and H > 0, fi > and H < 0, /i < and H > 0, and /i < 
and H < 0, respectively. In that Appendix we show that at fixed and finite values of the 
chemical potential and magnetic field the lowest energy Hamiltonian eigenstate of a family 
of states \r), S; r] z , S z ) with common values of r\ and S but different values of r\ z and S z is 
the [LWS,LWS], [LWS,HWS], [HWS,LWS], or [HWS,HWS] corresponding to the (I, I') sector 
choosed by the signs of the chemical potential and magnetic field. We have also calculated 
the smallest energy gaps relative to that state, which are given in Eq. (Cll). 

The main goal of this section is, however, to show that within all states with different rj 
and S values but the same eigenvalues r\ z and S z , i.e. of states belonging the same canonical 
ensemble, the ground state (18) has minimal energy. 
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Following the results of Sec. Ill, within all states I with the same eigenvalues r] z and S z 
the state of minimal energy has the form (18). Therefore, we can restrict our considerations 
to the set of non-LWS's and non-HWS's, \r), S; r] z , S z ), whose starting LWS's and (or) HWS's, 
\r), S; ±i], ±S), are "ground states" of the form (18). Other non-LWS's and non-HWS's with 
the same r] z and S z eigenvalues are of higher energy. 

The non-LWS's and non-HWS's belong the same canonical ensemble and, therefore, have 
common r\ z and S z eigenvalues. We emphasize that their starting states, \r), S; ±i], ±S), do 
not belong to that canonical ensemble. In the corresponding (1,1') = (sgn(r] z )l, sgn(S z )l) 
sector of symmetry U(l) ® U(l), the energy of the ground state (18) of eigenvalues r\ z and 
S z can be written as 

E° = E° so(4) (\ Vz \, \S Z \) + 2p Vz + 2fi HS z , (46) 

where E° so u^(ri, S) is the corresponding eigenenergy relative to the S'0(4)-Hamiltonian (2). 
Note that Eg , 4 Ji], S) is nothing but the term Eg ^ of the ground-state energy (Bll). 
Its i] and S dependence can be obtained from that expression by replacing the density and 
spin-density dependences by r\ and S dependences, respectively. The energy Eg , 4 Ji], S) 
reads 

N r™ 
E°so { A)(v,S) = ^j_Jk27rp° cm (k) [e(Q,,,, - \k\){-2tcosk - U/2} + U/4] , (47) 

where the function p Q c nm(k) is associated with the function p%^(v) through the integral 
equations (A33) and (A34) for the particular case of the ground-state distributions (B13) 
with Qiy defined in Eq. (A22). These two coupled integral equations have a unique solution 
which defines the functions p° c n^\(k) and p Q s n^i\(v). For fixed U the dependence of the energy 
(47) and functions p° c nii)(k) and p° s r U i\(v) on the r\ and S values is defined by the following 
normalization equations 

^ J_Jk2irp cm (k)[l - e(Q,,,, -\k\)]=v, (48) 

and 
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A/" r°° 

-^ J_Jv2<Kp° sm (v)[l - Q{B u ,/u -\v\)] = S. (49) 

These equations also define the dependence on r\ and S of the parameters Qiy and Biy of 
Eq. (A22). The use of Eqs. (13) and (14) reveals that (48) and (49) are equivalent to 
the normalization conditions (A36) and (A38) for the particular case of the ground-state 
distributions (B13). The above integral equations fully define the energy Eg Q uJr], S) and 
can be solved numerically. (A closed- form analytical solution is available for EgQt 4 \(Q, 0) 

10 

On the other hand, the energy E(r], S) of a non-LWS and non-HWS, \r), S; T] z , S z ), with 
the same values of r\ z and S z (but r\ z ^ ±r] and (or) S z ^ ±S) is, following Eq. (C2), given 
by 

E( V , S) = E° soi4) (rj, S) + 2 mz + 2fi HS z , (50) 

where the energy Eg uJr], S), Eq. (47), refers to the corresponding starting states 
\rj, S; ±r], ±S). Following our choice, these states are also ground states of form (18), but 
such that r] 7^ \r] z \ and (or ) S =fi \S Z \. 
We want to show that the energy gap 

E(rj, S)-E° = E 8O{4) (rj, S) - E° so(4) (\rj z \, \S Z \) , (51) 

where E°, E(r],S), and ^ 0(4) (r/, S) (and E° so{4) (\i] z \, \S Z \)) are given by Eqs. (46), (50), 
and (47), respectively, is positive for r\ z ^ ±r] and (or) S z ^ ±5. The energy Eg , 4 Jrj, S), 
Eq. (47), is a monotonous increasing function of both r\ and S, its minimum value being 
Eg O ,^(0,0). This refers to the "absolute" SO(4) ground state |]3|j20|]. Since 77 > \r] z \ and 



(or) S > \S Z \ [if S = \S Z \ (or r\ = \r) z \) we have that rj > \r) z \ (or S > \S Z \)], it follows that 
the gap (51) is always positive. Its smallest values correspond to the choices (a) rj = \r] z \ + l 
and S = \S Z \; and (b) r\= \r\ z \ and S — \S Z \ + 1. Let us evaluate the corresponding gaps 

A a = E° so{4) (\ Vz \ + 1, 15,1) - E° so{4) (\ Vz l \S Z \) , (52) 

and 
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A 6 = E° SO{4) (\ Vz \, \S Z \ + 1) - E° so{4) (\ Vz l \S Z \) , (53) 

respectively. The gaps (52) and (53) give the excitation energy of the states \\t) z \ + 
1, \S Z \; r/ z , S z ) and \\rj z \, \S Z \ + 1; r/ z , S z ), respectively, relative to the LWS's and (or) HWS's, 
I \Vz \, \S Z \; r] z , S z ), of the form (18). To evaluate these gaps we use the fact that in the ther- 
modynamic limit the excitation energy of the LWS's and (or) HWS's \\r] z \ + 1, \S z \;r] z + 
sgn(r] z )l, S z ) and \\r] z \, \S z \ + l;r) z , S z + sgn(S z )l) [also of form (18)] relative to the LWS and 
(or) HWS ||77 z |, \S z \;r) z , S z ) is zero. This excitation is the spin-flip ground-state - ground- 
state transition (i)-(iv) studied in Ref. [pOfl . Its excitation energy is of order l/N a and van- 
ishes in the present thermodynamic limit. (This is a condition for the continuous character 
of the magnetization curve defined by Eq. (42).) This implies that 

E°so(4)(\Vz\, \S Z \) + 2 mz + 2n Q HS z = £§ (4)(fel + 1, l&|) + 2/i[r/ 2 + sgn( Vz )l} + 2fi HS z , 

(54) 

and 

E°so(4)(\Vzl \S Z \) + 2fi Vz + 2fi HS z = E° so{4) (\r] z \, \S Z \ + 1) + 2 Wz + 2^E\S Z + sgn(S z )l] . 

(55) 

Taking into account the relation between the signs of /i (and H) and of rj z (and S z ) [see Eq. 
(29)], and combining Eqs. (52) and (53) with Eqs. (54) and (55), we finally arrive to: 

A a = 2|/i| ; A b = 2fx \H\ . (56) 

These are the smallest gaps of non-LWS's and non-HWS's belonging the same canonical 
ensemble (i.e. having the same eigenvalues i] z and S z ) relative to the corresponding ground 
state (18). Note that both (or one) of these gaps vanish(es) in the SO(4) sector [or SU(2) ® 
U(l) and U(l) ® SU(2) sectors] of parameter space. On the other hand, all Hamiltonian 
eigenstates with t] = and (or ) S — are both LWS's and HWS's of the corresponding 
algebras. This implies that in the canonical ensembles with r\ = and (or ) S = there are 
no non-LWS's and non-HWS's singlets of the eta-spin and (or ) spin algebras. 
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V. CONCLUDING REMARKS 



In this paper we have used the pseudoparticle-operator basis and perturbation theory- 



introduced in Refs. |L1,T3,m] to derive and study the ground states associated with all 
canonical ensembles belonging to the four U(l) (g> U(l) sectors of parameter space of the 
Hubbard chain in the presence of a magnetic field and chemical potential. Our results 
confirm the important role played by the pseudoparticle algebra in the low-energy physics 
of integrable quantum liquids: following the present study we find in Ref. [f2U| that the 
usual half-filling and zero-magnetic-field holons, antiholons, and spinons |3T| correspond to 
a limiting case of the general pseudoparticle representation. 

The simple form obtained for the ground-state, expression (18), has a deep physical 



meaning. It confirms fllTyi3|l and generalizes the fact that in the pseudoparticle basis the 
ground state of the many-electron quantum problem is a "non-interacting" pseudoparticle 
ground state of simple Slater-determinant form. This also holds true for canonical ensembles 



belonging to sectors of higher symmetry [20,28] and, therefore, in the pseudoparticle basis 



the ground state of canonical ensembles of all symmetries are always states I of that simple 
form. 

We have evaluated the energy gaps relative to the ground state of the non-LWS's and 
non-HWS's with common 7] z and S z eigenvalues. A complete proof of our ground-state 
expressions requires the calculation of the energy gaps of the states II with the same r\ z and 



S z eigenvalues [gS 



A more general Landau-liquid theory for the sectors U(l) ® U(l) including the states II 
can be constructed. These states can also be described in terms of pseudoparticles. However, 
in the sectors of lowest symmetry these requires, in addition to the pseudoparticles studied 
in this paper, new branches of "heavy" pseudoparticles |28 . 



The eight branches of a(l, I') pseudoparticles introduced in this paper have a deep physi- 
cal meaning. This is shown in Ref. |2(J where we relate the symmetry transformations of the 
set of pseudoparticles used in each canonical ensemble to construct the corresponding ground 
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state to the symmetry of that sector. In that reference we find that the pseudoparticles of 
this set always transform in the representation of the corresponding group of symmetry. 

Although the pseudoparticles associated with the states I are the transport carriers at 
low energy ]ID| and couple to external potentials |TI|, they refer to purely non-dissipative 
excitations, i.e. the Hamiltonian commutes with the currents in the subspace spanned by 
the states I [|l(| . Therefore, these pseudoparticle currents give rise only to the coherent part 
of the conductivity spectra, i.e. to the Drude peaks |]T0|JTl| , p!3[| . The finite- frequency part 
is associated with transitions involving the "heavy" pseudoparticles which are also needed, 
in the sectors of symmetry U(l) <8> U(l), to describe the states II f28 |. (In the Hilbert 



subspace spanned by those excitations, the Hamiltonian does not commute with the current 
operators.) 

As in the case of Landau's Fermi liquid theory pO| , |3Tf , the pseudoparticle perturbation 



theory uses as reference state the exact ground state of the quantum problem ||lT| , |T3f . Also in 
the construction of the above generalized Landau-liquid theory referring to the Hilbert space 
spanned by both the states I and states II [ ffijfl , the ground state which we have investigated 
and studied in the present paper plays a crucial role. 
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APPENDIX A: BETHE-ANSATZ SOLUTION OF THE FOUR 17(1) ® 17(1) 

SECTORS 

The BA solution associated with the Hamiltonian eigenstates I of the four (/, /') sectors 
is very similar to the solution of the (—1, —1) sector studied in Refs. [|3|j8|,|13| . 

In this Appendix we present the BA equations for the states I of the (I, I') sectors of 



symmetry U(l) (g> U(l) and introduce the two-pseudoparticle phase shifts [p~Q| , p~3] , [T~8|] and 
other quantities needed in the expressions presented in this paper. 

For each choice of N c q^ and N s (i t i>) occupied pseudomomenta values qj, Eq. (12), of 
the c(l, I') and s(l, I') pseudoparticles, respectively, (see the expressions of these numbers in 
terms of electronic numbers in Table 1) which describes one state I, Eq. (15), there is a set 
of iVc(J,J') re& l rapidity values, k- , and other N a {i^ real rapidity values, v ■ , which are the 



solution of the following N c ni>) + -^Ysftz') algebraic equations 

k 1 / =( lj + Jf Yj tan _1 (t$ - (l/u)smk 1 /) , j = 1, ..., N c{ i,i>) , (Al) 



"a At 



i'=i 



and 



9 N c(i,i') 

K 



*' =: "iw~ 53 tan 1 ( v J -(Vw)sinfc/ 
v a jf=1 



~ W "if tan- 1 ^ ( v »' _ v f)) , j = li ..., Nsm . (A2) 

o j'=l 

In the thermodynamic limit (A Q ,(; i //), iV a — >■ oo with n a (i,i') = N a (i t i/)/N a finite) the rapidity 
values kj and v- give rise to rapidity functions Kiy(q) and Sij>(q'), respectively, which 
are eigenvalues of the corresponding rapidity operators [|13|] . The set of algebraic equations 

(Al)-(A2) lead to the following two coupled integral equations 

1 ri {+) 
KiM =q + ~ AT'' dq'N sm (q') tan" 1 ^,^') - (V«) sin Jfy(g)) , (A3) 

and 



q s(l,l') 



l ^, (+) 



9 = - AT'' d l' N c{i,v){<l') tanrHSiAQ) - (V«) sin^i'^)) 

i /V +) i 

- - / : (M0 <tfN. m {ffi tan-^- (V(?) - V(</))) , (A4) 
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respectively, where the limits of the pseudo-Brillouin zones, q a ni>)i are given by Eqs. (19) — 
(21) with N a (i t ii) replaced by N^n^ and the pseudomomentum distributions, N c Q^(q), are 
the eigenvalues (and also expectation values) of the operators (26) relatively to the states I 
of form (15), ie 

N a (i,i'){q) = (i] z , S z \N a ( ltl ')(q)\7) z , S z ) . (A5) 

The pseudomomentum distribution N a ^i^(q) (A5) is 1 and for occupied and nonoccupied 
pseudomomenta, respectively, of the states (15). Therefore, the distributions (A5) fully 
define these Hamiltonian eigenstates I [see Eq. (27)]. We have that 

§ /5' ^ = N km > § /« l>) d ^m («) = *«<v) > ( A6 ) 

27Tj Ul>) 2 * J rJ,i>) 

where N*n t n and NaQ^ are the number of available «(/,/') pseudomomentum values , qj 
[see Eq. (12)], and numbers of a(l,l') pseudoparticles, respectively, given in Table 1. 

For each Hamiltonian eigenstate (15) there is one, and only one, pair of rapidity eigenval- 
ues Kiy (q) and Siy (q) . These are functionals of the pseudomomentum distributions. The 
solution of Eqs. (A3) and (A4) provides these rapidity functionals of the pseudomomentum 
distributions N a (i t i>)(q). 

It is easier to express the rapidity functions in terms of the eigenvalues of the normal- 
ordered operators (30) and (43), which define the pseudomomentum deviations. The rapidity 
functions can then be expanded in these deviations as jl3"Hl8| 

KiAl) = K${q) + K$(q) + K$(q) + ... (A7) 

and 

SiAti = S${q) + S$(q) + Sg!(q) + ... (A8) 

where K^Jlq) and SiJ{p) are the jth-order terms. Equations (A3) and (A4) allow the 
systematic evaluation order by order of all terms of the expansions (A7) and (A8). As shown 



in Ref. [13] for the case of the (—1,-1) sector, this deviation expansion corresponds to a 
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operator expansion in the pseudoparticle scattering order. The possibility of such expansion 



follows from the perturbative character of the pseudoparticle operator basis [|Tl],[T^,|T^] . 

Here we are interessed in the first-order terms of (A7) and (A8) which involve the two- 
pseudoparticle phase shifts pO|JT3| . [T8[1 . By using a recursion procedure, we find that the 
rapidities (A7) and (A8) may be simply written as 

K hV {q) = K^{QiAq)) , Sw{q) = S^tfVfa)) , (A9) 

where K\p{q) and S\jl(q) are the solutions that correspond to the choice of distribution 
(B3) of Appendix B (in this paper we want to confirm that this choice defines the ground 
state) and Qi,v(q) and Vi,i'(q) are functionals of the form 

QiAq) = q + Q?}{q) + Q?}{q) + -, (Aio) 

ViAq) = q + V§!(<1) + P$(?) + - , ( A11 ) 

which can be obtained by solving Eqs. (A3) and (A4) order by order. The results (A9)- 
(All) imply that the first-order terms of the rhs of Eqs. (A7) and (A8) may be written 



as 



^to) = ^f^cSh). (A12) 

and 

s m (q) = *J$M v m {qh (A13) 

respectively. We note that the functions dK( l ,{q)/dq and dS( J '(q) / 'dq obey the equations 

dK^iq) _ 1 

and 

dSJfi!{q) _ 1 



dq 2«fP am {S$hQJ) 
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(A15) 



respectively, where the functions 2Trp®n l ,\(k) and 2np s( , v Av) are the "ground-state" solutions 
of Eqs. (A33) and (A34) below. [The functions 2tx />°(_i -i)(^) an d ^np®, 1 _^(v) (with 
/°s(-i -l)^) = ucr ( u v), where v = A/it) are the usual ground-state distributions of Lieb and 
Wu @.] 

Solving Eqs. (A3) and (A4) to first order leads to 



a Jq a [l,V) 



2 (1) (?) = E AT" d ^ N M)^(Q, 80 > (A16) 



^ (1) (?) = e/ ( :, (,,,) w*fwW)*%{*,4) , (A17) 



a 9 a(U') 



where SN a ^ t i'\(q) = (r] z ,S z \ : N a (i t ii)(q) : |7fe, 5^). The four two-pseudoparticle phase shifts 
$^ tt , (g, </) can be written as 



•iTto,90 = «irf^^,^(rt|. (aw) 



< , fa8') = g , |gg- ) W, SiDj f' )(9 ' ) ), (A20) 



*!iW) = *!i' («$(«), S^tf))- (A21) 

Introducing the parameters 

Ql,V = K w \lFc(l,l')) > S, ( ,i/m = Sl°J{q Fsm ) , (A22) 

and 



.r 



_ sin Q u , 



yu> = B llV /u, (A23) 



(Q-1,-1 and -B-1,-1 are the usual cutoff parameters of the ground-state Lieb-Wu equations 
||) we find that the phase shifts $^ q/ obey the following integral equations: 
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^CC \ X 1 x 



I f y % d n gVj) 
Tri-yo, y l + (x-y") 



n\2 ' 



(A24) 



*^W) 



1 



TX 



tan (x — y ) + 



1 r». 



TX 



o y 1 + (x - y") 2 ' 



i,i' 



(A25) 



*£ (v,*') 



-— tan (y — x 

TX 



dy"G(y,y"M c (y",x' 



(A26) 



®/s (y, y') = - tan 



-i,V-V 



1_ r*l,, tajrW-i/) 



TX 



TX Z J-x° 



l + (y- x") 2 



+ r i 'Uy"G{y,y"W£{y" iy >) 



The kernel G(y,y') reads |TB| 



G(y,yO 



l 
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1 + ((y - y')/2)' 



2 V y - y 



where 



%) = - EO')taa _1 (y+KjO 



i=±i 



and 



1 



'(y) = - E0') ln d + (y + i4'> 

7T 



3=±1 



(A27) 



(A28) 



(A29) 



(A30) 



It is useful to introduce an alternative representation for the rapidity functions Kii'(q) 
and Sij>(q) of Eqs. (A3) and (A4) in terms of distributions p c {i,v){k) and p 8 (i,i')(v)- This leads 
to new equations which are equivalent to the latter equations. This second representation is 
less appropriate for the pseudoparticle operator basis but was, historically, the most widely 
used in BA problems PJ8|,[IS|]. The reason is that it leads to integral equations which, in 
some limits, are of mathematical standard type. However, and as we discuss below, the 
representation associated with the Eqs. (A3) and (A4) has a clearer physical connection to 
the BA operator algebra. 

Let us introduce the function p c (i,i')(k) such that 
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2np cm (K l Aq)) = l/[^^-]. (A31) 

This function is related to the distribution p s {i,i'){v) which is defined as 

2np sm (S l Aq)) = l/[^^-}. (A32) 



Equations (A31) and (A32) define the transformation q — > k for a(l, I') = c(l, I') and q — > v 

&) and K (<£$i'), 



for a(l,l') = s(l,l'), where the new variable k varies between K(q c J l ,- ) ) and K(q c J l ,- ) ) and 



the variable v runs from S(qJ ll n) to S(q s J l ^). 

Combining Eqs. (A3),(A4) and (A31),(A32) we find that p c n,i')(k) and p s (i,i>)(v) are the 
solutions of the following system of coupled integral equations 

IcosA; r s (9$ iir) ) ~ 27rp s(M / ) (i; / ) 

and 



vr u ysfaW ) ^ 'l + [* 



1 r K{ w jvat fiA 2 *pm)W 



-±r^j d y'N s{in (v>) 2 y )(t/ L ( A34 ) 

where N c Q t i^(k) and NgQ^^v) are the representation of the distributions iV c (^/)(g) and 
N s (i,i'){<l) (A5), respectively, in the fe,u space associated with the transformation g — > k 
[for a = c(7, 1')] and g — > v [for a = s(/, I')]. 

For (/,/') = (-1,-1) Eqs. (A33) and (A34) are similar to Eqs. (A4) and (A5), re- 
spectively, of Ref. |17| . However, Eqs. (A33) and (A34) are more general: the limits of 
integration of Eqs. (A4) and (A5) of Ref. |17j are only valid for Hamiltonian eigenstates 
differing by a small density of excited pseudoparticles from (18) (they are also valid for all 
eigenstates (15) where the occupied pseudomomenta of the functions (A5) are distributed 
symmetrically around the origin), whereas Eqs. (A33) and (A34) are valid in each canonical 
ensemble for all states I (15) in number of (17). 

For the particular case of a (—1, —1) state (18) the functions (A31) and (A32) are nothing 
but the usual distributions of Lieb and Wu PJT7fl. Equations (A33) and (A34) give the 
generalization for all states I (15). 
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Following Eqs. (A6), (A31), and (A32) the distributions p c (i,i')(k) and p s {i,v){v) obey the 
normaliztion conditions 

73-/ , ( _) , dk2irp cm (k) = N cm , (A35) 

Aa MUim) dk27C p cm (k)N cm (k) = N cm , (A36) 

^ / 5m ' ( "; ( /' )) dk2vp sm {k) = N: m , (A37) 



N a fSu'iijiW 



„(+) 

" dk2>irp sm (k)N sm (k) = N sm . (A38) 

[Here the numbers N^a^ and N a (i^ of the four (/, /') sectors of symmetry £/(l) ® C/(l) are 
given in Table 1.] 

As mentioned above, the representation associated with the distributions p c (_i i _i)(/c) 
and p s (-i-i) (v) was, until recently, the most used in BA in what concerns the description 
of states of form (18) and states whose distributions of BA quantum numbers differ from 
(18) by a vanishing density of these numbers PJ8|JT^]. One of the reasons for this is that 
from the mathematical point of view the integral equations (A33) and (A34) are, for these 
states, easier to handle than Eqs. (A3) and (A4). However, while the distributions (A5) 
of the latter equations are expectation values of the operators (26) and, therefore, have a 
clear physical meaning, we note that none of the distributions and functions p c n,i')(k) and 
p s (i,i'){v), N c (i,i'){k) and A s(M /)(u), and p c {i,v){k)N c {i,v){k) and p s (i,v){v)N s ^ v) {v) are expecta- 
tion values of any operator. (See, for instance, Fig. 3 of Ref. |32| where p c {-\ ,-i)(k) and the 



electronic momentum distribution are compared for zero magnetic field.) The latter distri- 
butions and functions are just a useful mathematical representation in the k, v space for the 
pseudomomentum distributions (A5) and eigenvalues Ki^(q) and Sij>(q) of the Hamiltonian 
eigenstates (15). 
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Finally, the expressions of the pseudoparticle bands (34) — (36) involve a third represen- 
tation for the two-pseudoparticle phase shifts (A18)-(A21) and (A24)-(A27) in terms of the 
k, v space variables. The corresponding phase shifts are given by 

5L//'/i ia J*//' [sink smk'\ , . nn . 

$'* (k, k') = &£ , , (A39) 

V u u 



^Uk,v') = ^(^y), (A40) 



* l £(v,h>) = * l £(v,*^p), (A41) 

and 

&/:(v,v') = ® l /:(v,v'), (A42) 

where the functions $^ Q / are defined by Eqs. (A24)-(A27). 

In Appendix B we study the energy of the Hamiltonian eigenstates I associated with the 
BA Eqs. (A3)-(A4) and (A33)-(A34). 
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APPENDIX B: GENERAL ENERGY EXPRESSION FOR THE STATES I 

In this Appendix we present and study the general energy expression for the states I with 
common eigenvalues i] z and S z corresponding to the (/, V) sectors of symmetry U(l) <8> U(l). 

The Hamiltonian eigenstates (15) are also eigenstates of the rapidity operator Ki^{q) 
and of the operators (26). Their energies can be written as 

E = E so{ 4) - HK m - N cm ] - iM)\H\[N; m - N, m ] , (Bl) 

where 



1 y a I 

v SO(4) 



Esn, • n = — /(! ( ) M ' ) dq [N cm (q){-2t cosK^(q) - U/2} + U/A] . (B2) 



27r J <^ 

Here N c ^^(q) is the c(l, I') pseudo momentum distribution (A5) and Ki >v (q) is the eigenvalue 
of the rapidity operator Ki t i'(q) defined by Eqs. (A3) and (A4): the rapidity eigenvalues of 
the states I of form (15), \rj z , S z ), are fully determined by the distributions (A5) through the 
BA equations (A3) and (A4). (In these equations Si t i>(q) is the eigenvalue of the rapidity 
operator Sij'(q).) 

For the case of the state (18) we denote the pseudomomentum distribution (26) by 

N a(i,i>)((i)- lt is g iven b y 

K(i,i>)(l) = (0,Vz,S z \N am (q)\0, Vz ,S z ) 

= ®(VFa(i,i>) - 9) , < q < q { +l n 

= ©(« - 9Fa(i,i0) ' «&) < « < ° " ^ B3 ) 

It is useful to express the energy (B2) in terms of the distribution p c (i,i')(k) of Eqs. (A33) 
and (A34). The term (B2) of the energy (Bl) can be rewritten in terms of that function as 
follows 

E so{4) = ^ f K ''" ^ dk27T Pcm (k) \N cm (k){-2tcosk- U/2} + U/A] , (B4) 

where N c (i,i')(k) is the distribution of Eqs. (A33) and (A34). The two energy expressions 
(Bl) with Eso(4) given by (B2) and (B4), respectively, are equivalent. For the sake of clarity 
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we have here expressed the term Eso(a) of the energy (Bl) in both the forms (B2) and (B4). 
These energies can be obtained by solving numerically Eqs. (A3), (A4), and (B2) or Eqs. 
(A33), (A34), and (B4). 

As mentioned in Sec. Ill, in each canonical ensemble of eigenvalues r) z and S z the state I 
of maximal energy corresponds to the pseudomomentum distribution (A5) of the particular 
form 

K(i,i')(Q) = (*,Vz, S z \N am (q)\*,r] z , S z ) 

= 1 - ©(<£(/,/') ~ 0) , 0<q< q§} tll) 

= 1 " 6(9 - q { -J m ) , q^ < q < , (B5) 

where 

a {±) -a {±) -a {±) (BG) 

Except for terms of 1/N a order, the pseudomomenta (B6) can be written as 

?£(M') ~ ±Q*a(i,i') , (B7) 

where 

71" 

q*a(i,i>) = q a (i,i>) - qFa(i,i>) = ^-[a* (/j/ ,) - a^mo] • ( B8 ) 

It follows that 

q*c(i,i') = Qc(i,i') — QFc(i,i') , q*s(i,i>) = q s (i,i') — Qfs(i,v) ■ (B9) 

The state I \*,r) z ,S z ) of the rhs of Eq. (B5) is of the form 

gM , q (+) t 

\*-,v z ,s z )= nt TT rf &Ui')]i*w>- ( m °) 

a=c,s (-) (+) 

q=q a( l,l>) 1=lJ(l,l>) 

Inserting the distributions (B3) and (B5) in Eqs. (A3) and (A4) defines the rapidities 
KiJ(q), Si l l(q) and K^*; (g), Sf], (g), respectively. (The rapidity functions K\j!(q) and S\i,{q) 
are also defined by Eq. (36).) We call the corresponding energies 
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Eo = ^o(4) - \»\lN* cm - N cm ] - pto\H\[N* sm - N. m ] , (Bll) 

and 

E. = E* so{4) - \fi\[N* cm - N cm ] - fi \H\[N: m - N sm ] , (B12) 

respectively, where E SO , A ^ and E* so , A ^ are given by (B2) with the c(/, /') pseudomomentum 
distribution defined by (B3) and (B5), respectively and the rapidities given by K u ,(q) and 
-K";;/(g), respectively. 

The energies Eg ^ and E* so , A ^ can also be written in the form (B4). Let us introduce 
the corresponding functions p%in(k), p^niiAv) and plai/\(k), p* s nirAv), respectively. These 
are solutions of the integral equations (A33) and (A34) with 

N° cm (k) = 9(Q M , - |*|) , N° m (v) = Q(B ltl ,/u - \v\) , (B13) 

and 

N* cm (k) = 1 - e(Q,V - 1*1) , N* sm (v) = 1 - e(B,V/ti - \v\) , (B14) 

respectively, where Qiy and B\y are defined in Eqs. (A22) and 

Qi,v = k $(Q*c(U')) > B tv/ U = S\*){q*,{iji>)) , (B15) 

respectively. In Eqs. (A22) and (B13)-(B15) we have used the values (22) and (B7)-(B9). 
For both the states (18) and (BIO) the limits of integration of Eqs. (A33) and (A34) are 
given by 

K$(±q cm ) = Ki;)(±q cm ) = ±vr , Si°J(±q sm ) = S$(±q. m ) = ±oc . (B16) 
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APPENDIX C: STATE OF MINIMAL ENERGY OF A FAMILY OF STATES 
WITH COMMON ETA SPIN AND SPIN 

In this Appendix we show that at fixed and finite values of the chemical potential and 
magnetic field the lowest-energy Hamiltonian eigenstate of a family of states \r), S; i] z , S z ) 
with common values rj and S but different eigenvalues r] z and S z is the [LWS,LWS], 
[LWS,HWS], [HWS,LWS], or [HWS,HWS] corresponding to the (/,/') sector choosed by 
the signs of the chemical potential and magnetic field. 

Obviously, we have that 

H S o(4)\v,S;t] z ,S z ) = E SO (4)(v, S)\r], S;r] z , S z ) , (CI) 

where Eso(4)(Vi S) * s the eigenenergy of the whole family of states of eta spin 77 and spin S 
corresponding to the Hamiltonian (2). Moreover, the eigenenergy of \r], S; i] z , S z ) relative to 
the Hamiltonian (1) is 

E( Vz , S z ) = E SO (4)(v, S) + 2fi Vz + 2fi HS z , (C2) 

whereas the eigenenergy of the four above LWS's and (or) HWS's is 

E[lws,lws] = E SO (4) (v, S) - 2/j,r] - 2/j, HS , (C3) 

E[lws,hws\ = Eso(4) (v, S) - 2/ir) + 2/i HS , (C4) 

E[hws,lws] = Eso(4) (v, S) + 2/ir] - 2fi HS , (C5) 

and 

E[hws,hws] = E SO (4) (v, S) + 2/iT] + 2/j, HS . (C6) 

For i] z 7^ ±77 and (or) S z 7^ ±S out of the four energies 

E( Vz , S z ) - E [LWS , LWS] = 2fi( Vz + v ) + 2fi H(S z + S) , (C7) 
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E( Vz , S z ) - E [lws ,hws] = MVz + V) + 2fi H(S z - S) , (C8) 

E( Vz , S z ) - E [HWS , LWS] = 2fi( Vz - rj) + 2pt H(S z + S) , (C9) 



and 



E(r) z , S z ) - E [hws> hws] = 2/i(^ - rj) + 2/i H(S z - S) , (CIO) 

the energy difference corresponding to the state \rj, S; —i], —S), \rj, S; —rj, S), \rj, S; rj, —S), or 
\t],S;t],S) which belongs to the (1,1') sector choosed by the signs of the chemical potential 
fi and magnetic field H is positive and maximal. This follows directly from the form of the 
four energies (C7)-(C10). 

For example, consider that \x > and H > 0. In this case the eta-spin and spin 
projections are such that i] z < and S z < [see Eq. (29)], respectively, and we are 
in the (—1,-1) sector of symmetry U(l) <E> U(l). It follows then from Eqs. (C7)-(C10) 
that the [LWS,LWS] Hamiltonian eigenstate \rj, S; —77, — S) has minimal energy and the gap 
E(rj z ,S z ) — E[ LWS ,lws] is positive when rj z 7^ —77 and (or) S z 7^ —S. (Also the energies 
E[lws,hws\ — E[lws,lws], E[ HWStLWS ] — E[ LW s,lws\, an d Ey HWS) Hws\ — E[ LW s,lws] are in 
this case positive.) There are two minimal gaps, A,, and As, which correspond to eta spin 
and spin, respectively. They refer to the states 1 77, S; —i] + 1, —S) and \r], S; —rj, —S + 1), 
respectively. Following Eq. (C7) they read A^ = 2/i and A5 = 2fi Q H, respectively. 

The same applies for other signs of the chemical potential and magnetic field, i.e. for the 
remaining three U(l) ® U(l) sectors. The corresponding LWS and (or) HWS is the state of 
minimal energy of the family of states \i], S; i] z , S z ). In the general case, the minimal gaps 
are given by 

A„ = 2H; A s = 2 f i \H\. (Cll) 
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(-1,1) (1,-1) (1,1) 

N 2N a -N 2N a - N 

N a N a N a 

iV T N a - 7V T N a - 7V| 

iV| N a -N t N a - iV T 

2k F 2[7T-k F } 2[TT-k F ] 

71 71 7T 

k F i 7r — k F ^ 7X — kpi 

k F { vr — kpi 7T — kp^ 



TABLE 1 - Values of the numbers N a ^ t i^ and N*n m and of the pseudo- Fermi points and 
limits of the pseudo-Brillouin zones (22) in the four (/, /') sectors of symmetry U(l) <g> U(l). 





I- 1 


N c (i,i') 


N 


N kn 


N a 


N s (i,i>) 


Ni 


N* 
s(M') 


iV T 


QFc{l,l') 


2k F 


Qc(i,i') 


71 


QFs{l,l') 


kFi 


Qs(i,i') 


kp^ 
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